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We have measured the complex surface impedance of 2H-NbSe2 in the mixed state over a wide 
range of magnetic field(0-2T) and frequency(10-3000MHz). A crossover between pinned and viscous 
dynamics of flux-lines is observed at a pinning frequency, u p . The measured ui p is compared to that 
predicted by a single-particle "washboard potential" model for the pinning interaction. When the 
flux lattice is prepared in an ordered state by zero-field cooling or by driving it with a large current, 
uip is in good agreement with the predicted value. However, when it is prepared in a metastable 
disordered state by field cooling, uj p is nearly two orders of magnitude higher than expected indicating 
a complete break-down of the single-particle model. 



PACS numbers: 74.60.Ge 74.60.Jg 74.60Ec 

Pinning of the magnetic flux line lattice (FLL) due to 
material disorder plays an important role in the trans- 
port properties of type II superconductors |Q. It can 
lead to non-dissipative DC transport and to enhanced 
screening of electromagnetic fields. Most studies of pin- 
ning have focused on the DC critical current, j c , at which 
the FLL breaks loose from the pinning centers and starts 
moving. However, the observation of phenomena such 
as current induced annealing and finite response 

times to large current pulses pj, have shown that the 
threshold measurements do not completely characterize 
the pinned state. A complementary approach is to probe 
the frequency dependence of the response to subcritical 
AC currents which induce small oscillations about the 
pinned state. This is described in a mean-field model 
by the equation of motion: 



rju + K p u = $oj{t) 



(1) 



where j(t) — je~ luJt is the applied current density (in 
general, thermal driving forces can be included), u is the 
displacement from equilibrium, n p is a restoring force 
constant due to pinning, 77 is the viscosity and <E>o = ff 
the flux quantum. At a characteristic pinning frequency, 
ujp = K„/rj there is a crossover from a low frequency 
regime where the FLL is pinned and the response non- 
dissipative to a high frequency regime of free flux motion 
with viscous response. Such a crossover was indeed ob- 
served in early AC resistivity measurements of alloy su- 
perconductors P], but recent results in high T c samples 
were not compatible with mean-field behavior S. 

In this Letter we report on swept frequency measure- 
ments H (10-3000MHz) of the surface impedance in the 
mixed state of the low T c superconductor 2H-NbSe2. The 
results are in good agreement with mean-field calcula- 
tions of the electromagnetic response Q| which include 
contributions from Cooper pairs, the normal fluid, flux 
flow and pinning. A distinct crossover frequency is ob- 
served separating the regimes of pinned and viscous re- 
sponse. We find that u> p exhibits a striking sensitivity to 



the state of the FLL: for the same field and temperature 
the FLL can have pinning frequencies which differ by as 
much as two orders of magnitude, depending on how the 
FLL is prepared. This surprising result is interpreted in 
light of pronounced metastability effects that were ob- 
served in DC and pulsed current measurements [^[^J. If 
the FLL is prepared by zero field cooling to T < T c and 
then applying a field, H, the FLL enters one of two sta- 
ble states. For (H,T) below a well-defined transition line, 
T m (H), which is marked by a jump in j c , the state is or- 
dered and weakly pinned, while above T m (H) it is more 
strongly pinned and disordered. However, cooling the 
FLL from above T c in a constant field always results in 
a disordered FLL. Below T m (H), the disordered FLL is 
metastable: it can be annealed into the ordered state by 
driving it with a current that exceeds j c , by changing the 
field or by mechanical shock (after annealing j c is up to 
six times lower). The ordered state is stable against such 
variations. Our present results show the dramatic differ- 
ences in the AC dynamics of the two states below T m (H): 
in the ordered state the FLL has a relatively low lu p which 
is consistent with the measured j c in a single-particle 
washboard model. By contrast, in the metastable disor- 
dered state, ujp is nearly two orders of magnitude higher 
than predicted by this model, indicating the existence of 
complex underlying dynamics. 

The experimentally accessible quantity at high fre- 
quencies (where the AC fields are shielded out of bulk 
samples) is the surface impedance, Z s = R s — iX s |jTo| . 
R s and X s are the surface resistance and reactance. In 
the mixed state, Z s can be expressed in terms of a com- 
plex penetration depth, A, by Z s = —i^uoK. A is deter- 
mined by the response of Cooper pairs, the normal fluid 
and flux- lines [Q: 
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Here A is the London penetration depth, 5 n f is the 
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normal fluid skin depth, S v = (2p v / pqlo) 1 ^ 2 is the 
flux-line skin depth and p v is the flux-line resis- 
tivity. Our results are expressed in terms of the 
reduced variables, r s = R S (T)/R S (T C ) = Im(§^) and 
x s = X S (T)/X S (T C ) = Re(j±), which eliminate a trivial 

lo 1 / 2 dependence, due to the skin effect. The normal state 
skin depth, 6 n , is almost independent of H and T in our 
samples. This was determined from measurements of the 
normal state resistivity, p n . Most results described here 
are in a range of H,T, and 10 where A <C 5 V , S n f and the 
response is dominated by the flux- lines: 



r s - ix s = (-2ip v / p n ) 1/2 . 



(3) 



If we assume the Bardeen-Stephen value for the viscosity, 
?; = &qH C 2/ p n , with H C 2 the upper critical field, then in 
the model described by ([!]), the flux-line resistivity is: 



p v {u) = p v i(u>) - ip V 2(u) 



iuiujp H 



oj 2 +oj 2 H, 
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(4) 



In the free flux flow limit (uj 3> lu p ), the response is fur- 
ther simplified and reduces to that of a normal metal 
with a purely real resistivity, p n H / H C 2(T): 

r s (T,H) = x s (T,H)=(H/H c2 (T))? T<T C (H). (5) 

Measurements were done on two single crystal sam- 
ples. For sample 1 (~2.5 x 2.5 x .1mm) T c = 7.2K, 
AT C =80mK, and the residual resistivity ratio, RRR==23. 
For sample 2 (-4 x 4 x .025mm) T c = 5.85K, AT c =80mK 
and RRR =9. The parameters indicate that sample 2 is 
relatively impure, and that both have good homogene- 
ity. The field was along the sample's c-axis, and the AC 
current was in the a-b plane. We will focus on data for 
sample 2. 

The experimental setup J^] is shown in Fig. [l]. The 
sample is mounted on a sapphire plate, which is attached 
to a stainless steel holder inside a vacuum can. A 50S1 
coplanar transmission line consisting of three Cu strips 
(one central conductor between two ground planes) on 
an alumina substrate is placed close to the sample. The 
sample's surface resistance (reactance) introduces a slight 
change (< 10~ 4 ) in the magnitude (phase) of the sig- 
nal propagating along the transmission line. This ef- 
fect is isolated from a large background by modulating 
the sample temperature at a low frequency (~5-10Hz) 
with a resistive heater. The mixer puts out a signal 
at the modulation frequency, which is proportional to 
cos((f>)dR s /dT + sin((j))dX s /dT, and is detected with a 
lock-in amplifier. By tuning 0, the phase difference be- 
tween the LO and RF ports, the derivative with respect 
to temperature of R s and X s can be measured separately. 
All the data were taken in the linear regime: the RF cur- 
rents in the sample were ~lmA (j ~ lA/cm 2 <C j c ) and 
the amplitude of the flux-line motion (<.lA at 10MHz) 
was negligible compared to their spacing. 
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FIG. 1. Upper Panel: schematic diagram of the measure- 
ment technique. Lower panel: cross sectional view of the 
sample holder and transmission line. 



We first present the results in the disordered state. The 
FLL was prepared in this state by cooling from above T c 
to 3K in a constant field. In Fig. || we show the T de- 
pendence of r s and x s . The data was obtained by slowly 
ramping T as dR s /dT or dX s /dT was measured. r s was 
obtained from the raw data in the following manner. The 
zero field data were integrated, with the integration con- 
stant fixed by assuming that R s =0 for T <C T c . The inte- 
grated data were scaled to be equal to 1 at T c . The value 
of the scaling factor depends on the coupling between 
sample and transmission line which is independent of H 
and therefore the same scaling factor was used for the fi- 
nite field data. The finite field integration constants were 
set by using the fact that Z s is independent of field in 
the normal state. The same basic procedure was used for 
the reactance data, but determining the integration con- 
stants was complicated by the fact that even at low tem- 
peratures and zero field, X s has a finite value p,QUi\ and 
the absolute value of A was not measured directly. The 
necessary parameter xo = x s (3K, OT) = 2\(3K,0T)/5 n , 
was determined by the method described below. 

In Fig. |a the measured r s (T) at 2.8 GHz is com- 
pared to the calculated values. For fields >.1T, the 
approximate expression (|^) is in good agreement with 
the data. This calculation had no adjustable parame- 
ters since H C 2(T) was obtained from DC measurements 
of T C (H). The agreement indicates that for this fre- 
quency the FLL is in the free flux flow regime and con- 
sequently the pinning frequency for the disordered state 
LUp -C2.8GHZ. At lower fields, S v ~ A and the approxi- 
mate expression is no longer adequate. Fitting the data 
at .01, .02, and .05T to the full expression (||), with xq as 
an adjustable parameter, we obtain xo ~.09 at 2.8GHz. 
This value was used to process the reactance data shown 
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FIG. 2. The temperature dependence of the surface resis- 
tance and surface reactance at 2.8GHz (a and b) and 13MHz 
(c and d) for field (□), .02T (O), -IT (•), .25T (A), .5T 
(V), IT (O), and 1.5T (<g>). The FLL was field cooled in all 
cases. The dotted lines were calculated using (^|); the solid 
lines were calculated using (^). The double arrows show the 
deviation from free flux flow behavior. 

in Fig. |^b. We see that the full expression is in good 
agreement with the data at all fields for both r s and x s . 
(A small contribution from pinning, determined by the 
measured pinning frequencies, was included in these cal- 
culations. At the higher fields, this is the main source 
of the slight difference between the calculations using (|^) 
and (||)). Since xq = 2X ( 3 ^,ot) ^ depends only on 
p n , which was measured, the value of xq can be used 
to determine A. From the data at several frequencies in 
the GHz range, we estimate A ~1200-1400A in sample 
2 and ^1000-I200A in sample 1. Previous estimates are 
scattered over the range 700-2500A 0. 

In Figs. ||c and ||d we show the results at 13MHz. In 
this case, both r s and x s are well below the free flux flow 
values, indicating that cup 1 > 13MHz. The frequency de- 
pendence of the response over the entire range, shown in 
Fig. |^, exhibits a crossover from pinned to free behavior. 
The best fit of the data to (0), with p v given by (||), is 
obtained for ujp = 125MHz. The upper inset of Fig. ||a, 
shows the temperature dependence of tup 1 and j c at .5 T. 
The lower insets in Fig. 3 show that the data for all fields 
collapse when plotted as a function of uj/u)p(T, H). This 
universal behavior demonstrates that in this state the 
frequency dependence of the response of the FLL is com- 
pletely characterized by one parameter, cup 1 , as expected 
from (1). 

We now turn to the response in the ordered state. 




Frequency MHz 

FIG. 3. The frequency dependence of R a (a) and X 3 (b) at 
3K and 0.5T. (a), upper inset: field and temperature depen- 
dence of the pinning frequency. Lower insets: the collapse of 
the data for all fields. Symbols are the same as in Fig. g. 

When the FLL is prepared in this state and then heated 
above T m , it undergoes an order to disorder transition 
that is irreversible in temperature Thus, dR s /dT and 
dX s /dT cannot be obtained with this technique over the 
entire temperature range and it is not possible to set 
the integration constants (since this requires data up to 
T c ). Although the absolute values of r s and x s cannot be 
determined, the relative changes, Ar s (T) and Ax s (T), 
are still accessible for T < T m . At frequencies in the 
GHz range, the results are independent of how the FLL 
is prepared. In Fig. |i]a, we plot Ar s (T) at 13MHz in 
both states. The data in the disordered state is essen- 
tially constant in temperature, consistent with the high 
pinning frequency of this state. In the ordered state how- 
ever the data has a substantial slope which is close to its 
value at 2.8GHz, indicating that the pinning frequency 
in this state w° <13MHz. 

Since the value of lo° appears to be in the MHz range, it 
can be measured by more conventional means. At these 
frequencies the sample thickness is much less than 5 V , so 
that the skin effect is negligible and p v could be accessed 
directly with a 4-lead technique. The sample used for 
these measurements was from the same batch as sample 
2 and had nearly identical parameters. A background 
signal was determined from zero field measurements at 
T <C T c and subtracted from the data. In Fig. ^b we see 
that in the disordered state both p v \ and p v i are van- 
ishingly small over the entire frequency range, 0-200kHz, 
indicating that ujp 3>200kHz, in accord with the previous 
result. By contrast, p v \ and p V 2 are finite in the ordered 
state and their frequency dependence gives , uj° ~1MHz. 

To interpret the results, we use a model where the 
FLL is treated as a single particle moving in a washboard 
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FIG. 4. Comparisons of the AC response for the disordered 
and ordered states below T m (a) Ar s (T) at 13MHz. (b) The 
frequency dependence of the complex resistivity. The noise 
near 75KHz is instrumental. 



potential: V(u) — Vq(1— cos(^ i )). u is the displacement 
from equilibrium and r p is a characteristic length scale for 
the interaction of the FLL with the pinning sites. In this 
model, the critical current is j c — and the restoring 



force constant is k z 



V ( 



2tt\2 
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0*0 
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was previously observed with DC f|,[|[L2| and pulsed cur- 
rent measurements (J] , as well as in computer simulations 
]l3[ . However, in the AC measurements, no rearrange- 
ments can take place because of the low driving currents 
and the short time scales. 

In summary, the AC response in both states exhibits a 
crossover between pinned and free behavior which can be 
described by a simple equation of motion. In the ordered 
state the connection between lo p and j c is accurately pre- 
dicted by a single particle washboard model, but for the 
metastable disordered state the model breaks down. The 
two states discussed here are possible candidates for the 
proposed vortex glass (disordered) and pinned lattice or 
Bragg glass (ordered) Q . Our results indicate that each 
of these states can be supercooled or superheated across 
T m into the stability region of the other In this 

picture, the metastable disordered state corresponds to 
a supercooled vortex glass. 

We thank N. Andrei, D.A. Huse and R. Walstedt for 
useful comments. Supported by NSF-DMR-9401561. 



leading to: 
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2tt j c p„ 



(6) 



For the ordered state data shown in Fig. Wb, j c = 
70A/cm 2 . Choosing r p —ao (where ao = (Qo/Hy-/ 2 is the 
flux line spacing), as is usually done, gives cj°=1.4 MHz, 
in rough agreement with the data. In the disordered 
state (^(3^,0.57) = 350A/cm 2 ), the washboard model 
gives uj p = 2.6 MHz, which is 45 times lower than the 
observed value. This discrepancy is similar everywhere 
in the H-T plane except above T m (H) (where the disor- 
dered state is stable), in which case it is much smaller. 
The washboard model can give a higher oj p if r p is re- 
duced. But even using the smallest length scale in the 
problem, the coherence length, and taking into account 
the structure of the flux-line core [Q leads to values that 
are not much larger: w p (3K,.5 T)=4.7 MHz. In addition, 
the washboard model cannot account for the qualitative 
difference in the temperature dependences of lj p and j c . 
Thus, this model successfully describes the connection 
between j c and lo p in the ordered state, but fails to do 
so in the metastable disordered state. This is due to the 
single particle treatment of the FLL in which both j c 
and ujp are determined by the shape of a pinning poten- 
tial. In this model, j c is the current to drive the entire 
FLL out of a potential. However, the metastable state 
reorders as it depins and j c for this state may actually be 
the current at which the FLL becomes unstable to rear- 
rangements which create continuous channels of mobile 
flux lines. Evidence for channel formation and growth 
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